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I. INTRODUCTION 

The resistivity of a normal metal (N) on distance scales 
greater than the mean free path £ is characterized by 
Drude theory. Because the Drude model is based on the 
diffusive motion of electrons, the resistance scales with 
material or probe geometry. On length scales shorter 
than £, electrons move ballistically and the resistance 
is not so simply related to the geometry. This ballis- 
tic regime, which is becoming increasingly relevant with 
the miniaturization of electronic devices, can be accessed 
by tunneling into a material from a contact of size a <C £■ 
In such contacts, rich spectroscopic information may be 
obtained from measurements of differential conductance 
(d//dU) or its derivative (d 2 I/dV 2 ) with respect to ap- 
plied voltage. For example, in a normal-normal junc- 
tion (N/N), a measurement of the conductance derivative 
yields information about the spectral weight of interac- 
tions, which are typically electron-phonon in nature. 1 ' 2 

New phenomena occur when the point-contact is be- 
tween a normal metal and a superconductor (N/S). For 
example, if there is no barrier between the two metals, 
the current through the junction increases by a factor of 
two below a certain energy A. This energy is associated 
with the formation of Cooper Pairs in the superconduc- 
tor. On the other hand, for a large barrier the current 
through the junction goes to zero for energies less than 
A. 

Point-contact spectroscopy of N/S junctions has been 
predominantly used to measure the energy gap and as- 
sociated density of states. Using this information, lower 
bounds on the Fermi velocity and coherence length (the 
spatial size of a Cooper pair) of the superconductor can 
be obtained. Measurements on superconductors to ob- 
tain such information have been performed on materials 
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FIG. 1: (Color online) Schematic of a point contact measure- 
ment where a normal metal tip touches a superconducting 
sample. A differential conductance measurement is made by 
applying a stepped DC current with small AC modulation 
across the junction while the AC voltage is measured. 



such as niobium (Nb), 3 ~ 5 cuprates, 6 Sr 2 Ru04, 7 MgB 2 , 8 
and heavy fermion materials. 9 In this paper, we present 
an experiment, suitable for an advanced undergraduate 
laboratory, to study the N/S spectroscopic behavior of a 
point-contact between a gold (Au) tip and Nb film (see 
Fig. 1). 

In order to understand how reflection and transmission 
occur at the interface between a normal metal and a su- 
perconductor, it is important to understand some basic 
facts about superconductivity. At the transition temper- 
ature of a superconductor T c , the resistance drops to zero 
and an energy gap A ~ 1CU 3 eV, opens up in the den- 
sity of states. As explained by Bardeen-Coopcr-Schricffcr 
(BCS) theory, 12 this energy gap — also known as the pair- 
potential energy or binding energy — arises because it is 
thermodynamically more favorable for electrons of oppo- 
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site spin to bind together into a spin-zero boson called a 
"Cooper pair." Hence at T — 0, when there are no ther- 
mal fluctuations, the electrons will exist only in Cooper 
pairs. 

In N/S point-contact transport measurements, three 
major processes occur at the interface: (1) transmission, 
(2) reflection, and (3) Andreev reflection. All can occur 
to varying degrees, although differences in the size of the 
junction and the tunneling barrier can alter the prob- 
abilities for any given event to occur. The "height" of 
the barrier at the N/S interface can be characterized by 
a dimensionless parameter Z; a low barrier that is easy 
to cross corresponds to Z sa while a high barrier that 
inhibits current is characterized by Z > 1. The resulting 
differential conductance Gns for high and low barriers at 
T = is shown in Fig. 2. 

Let us now return to a description of the three pro- 
cesses mentioned above. For energies above the gap 
[E > A), Cooper pairs break and form single parti- 
cle states in the superconductor. Such states are called 
"quasiparticles" because a particle and its influence on 
the local environment move together as a single particle- 
like entity. In this energy regime, process (1) transmis- 
sion occurs. Here, single electrons are transmitted into 
the superconductor and manifests itself as a linear I-V 
curve (constant differential conductance). 

For energies below the gap (E < A), the differential 
conductance depends on the constriction and tunneling 
barrier. In the case of a large constriction where trans- 
port is not in the ballistic regime, process (2) reflection 
occurs independent of whether the tunneling barrier Z 
is small or large. In this case a single particle fermion 
cannot be transmitted into the superconductor which is 
composed of two-particle bosonic states. The electrons 
from the normal metal therefore reflect from the interface 
of the superconductor and hence the current and dl/dV 
are zero for this energy regime. Right at the gap energy 
A, there is a discontinuity in the I-V curve that gives rise 
to a peak in dl/dV. The width of this peak is tempera- 
ture dependent as broadening can occur due to thermal 
fluctuations that destroy the Cooper pairs. 

If the constriction is small enough to be in the ballistic 
regime, then for a small tunneling barrier with E < A, 
process (3) Andreev reflection occurs and the differential 
conductance increases by a factor of two from the conduc- 
tance above the gap (see Fig. 2a). Physically, this process 
can be characterized as shown in Fig. 3, where in order 
for a spin-up electron that is injected with energy E < A 
to be trasmitted, it must form a Cooper pair with an 
electron of opposite spin in the superconductor. Conse- 
quently, a hole — a mathematical construction to describe 
the absence of an electron — is reflected back with oppo- 
site spin and momentum of the incident electron. Hence 
if the incident electron carries a current evp, where v-p is 
the Fermi velocity, then a current of — ev-p is reflected as 
a hole. The net result is a current of 2et>F in the normal 
metal and a current of 2eup in the superconductor as a 
result of the Cooper pair. 
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FIG. 2: (Color online) Normalized differential conductance 
at positive bias voltage in the (a) Andreev and (b) tunneling 
regimes. In Andreev reflection the differential conductance 
increases by a factor of two in comparison to the normal state 
above the gap energy A. 
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FIG. 3: (Color online) A diagram showing Andreev reflection, 
where an electron in the normal metal with energy E = Ef + 
eV (and eV < A) is incident upon a superconductor. Note 
that in the body of the paper, Ef is referred to as E — 
because only excitations above this energy are considered in 
the theory. 



In a real N/S contact, a tunneling barrier will exist 
and the resultant differential conductance measurement 
will be somewhere in-between the two extreme tunnel- 
ing barrier cases of Fig. 2. However, in the experiment 
presented in this paper, by using a differential screw to 
gradually press an Au tip into superconducting Nb, the 
tunneling barrier Z can be varied over a range of values. 

This paper is organized as follows: In section II we 
discuss the Blonder-Tinkham-Klapwijk (BTK) theory 3 
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FIG. 4: (Color online) Potential energy diagram for the BTK 
model. There is a constant potential energy A in the super- 
conductor and a barrier at the N/S interface modeled by a 
5-function of strength H. 



and derive the differential conductance and other rele- 
vant physical quantities. In section III we discuss our 
differential screw-dipping probe apparatus and the cir- 
cuit used to obtain the conductance curve as a function 
of the bias voltage applied across the junction. Finally, 
in section IV we present the results of our own differen- 
tial conductance measurements, including an estimate of 
the energy gap and lower bounds for the BCS coherence 
length and Fermi velocity of Nb. 



II. THE BTK MODEL 

The BTK theory used to model the point-contact junc- 
tion considers the transmission, reflection, and Andreev 
reflection of a lone electron at the boundary between a 
normal metal and a superconductor. 3 We take the po- 
tential to be a (^-function barrier of amplitude H at the 
interface, and a constant A inside the superconductor 
(see Fig. 4). The incident wave is assumed to be a one- 
dimensional (ID) plane wave traveling from the normal 
metal into the superconductor. A ID model is sufficient 
for Nb as its gap is known to be spherically symmetric 
in momentum space. We first calculate the probabilities 
for reflection and transmission of an electron at the N/S 
interface using familiar techniques learned in quantum 
mechanics. Afterwards, we use the results to determine 
an expression for the macroscopic current that will be 
observed in our apparatus. 



A. Solution of the BTK Model 

In order to find the probabilities of the various 
processes, we solve the Bogoliubov-de Gennes (BdG) 
equations, 3 which are slightly modified versions of 
Schrodinger's equation. If we denote the wave function 



f(x,t) 
g(x,t) 



(1) 



where f(x,t) and g(x,t) represent the electron and hole 
wave functions respectively, then the BdG equations are 



2m dx 2 
A(x)a x tp. 



fi(x) + V(x) ) a z tp 
(2) 



Here, a x and a z are Pauli spin matrices 
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0(x) is the step function 
Q(x) = 



x>0 
x < 



(3) 
(4) 

(5) 



V(x) is the potential energy, (j,(x) is the chemical poten- 
tial, and A(x) represents the gap. The rate at which 
quasiparticles decay to the BCS ground state inside the 
superconductor is represented by the parameter T; if t is 
the lifetime of these particles, then 



h 

T 



(6) 



The effective Schrodinger-like equation for this system 
looks like a generalized version of the traditional prob- 
lem of quantum tunneling through a ^-function barrier 
with a few extra terms. The parameter T governs quasi- 
particle decay through scattering and/or recombination 
effects. 14,15 In practice, T serves as a phenomenological 
parameter that broadens the peak shown in Fig. 2b. The 
chemical potential is the energy required to add a par- 
ticle to the system; since both electrons and holes are 
fermions, this is equivalent to the Fermi energy of the 
system at T = 0. The presence of the Pauli matrices is 
a convenient way to represent electrons, holes, and their 
interactions into one equation. Phenomenologically, the 
term Aa x allows for a coupling between electrons and 
holes. 

In the BTK model, we consider the case where fi(x) — 
[i and A(x) = AQ(x) for constant fi and A, and V(x) — 
H8{x). To characterize the height of the potential barrier 
at the interface, it is convenient to define a dimensionless 
parameter 
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where i>p s is the Fermi velocity of the superconductor. 
To solve for the steady-state plane-wave solutions at the 
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N/S interface, we use trial wave functions of the form 
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These incident, reflected, and transmitted waves lead to a 
more complex version of a common undergraduate prob- 
lem on transmittance through a <5-function barrier. 

There are several important points to make regarding 
the form of tp. Firstly, there are four distinct wave num- 
bers k and q . The reason to distinguish between fc's 
and g's is because the inclusion of Q(x) alters the disper- 
sion relation u>(k) inside of the superconductor. Also, the 
signs on the diagonal of a z require that we distinguish be- 
tween the + and — wave numbers. Secondly, inside the 
superconductor the waves are no longer pure electrons 
and holes but instead a combination of both — they are 
quasiparticles. Because the Hamiltonian is no longer a 
diagonal matrix in the electron-hole basis, its eigenvec- 
tors will have two non-zero components. For any given 
energy there are two quasiparticle states and we denote 
the one with a positive (negative) group velocity by k + 
(k~). Finally, a reflected hole propagates with a wave 
number of +q~ and not — q~ . Because a hole is the ab- 
sence of a particle, differentiating the dispersion relation 
gives a group velocity with the opposite sign. One of the 
quasiparticles will propagate forward with a wave num- 
ber of — fc~ for the exact same reason (see Fig. 11.7 of 
Ref. 12 for a nice picture of this dispersion relation). 

The boundary conditions for this problem are continu- 
ity of the wave function 



^n(0)=Vs(0), 



(11) 



where -0s = V't and t/>n = tpi + , and a discontinuity in 
the slope of the wave function across the <5-function 
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Using these conditions, we can derive the relations 

hq ± = y/2m(n±E), (13) 



Hk* = J 2m n± J(E + iTf - A 



and 
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We can solve for a and 6 from Eq. (9) if we approx- 
imate q + = q~ = k + = kT = fcp (or equivalently 



/i E, A, r). Because we are dealing with energies 
of only a few meV whereas a typical (normal) metal has 
a Fermi energy of 1-10 eV, this is a reasonable approxi- 
mation for any normal metal (see Table 2.1 of Ref. 10). 
We then obtain 
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where 



6 = - 
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(16) 
(17) 

(18) 



Using these results, we can calculate the probabilities 
of Andreev reflection A(E) = aa* and normal reflection 
B(E) = 66*. Lastly, because Nb and Au have different 
Fermi velocities, we must use an effective Z e g that satis- 
fies 



Z 



eft 



Z A 



4r 



(19) 



where r — vf n /vf s is the ratio of the Fermi velocities in 
the normal and superconducting metals. 16 By using Z e g 
we are able to place a lower bound on the Fermi velocity 
of the superconductor. 



B. Differential Conductance 

We now derive the differential conductance of a N/S 
tunnel junction from first principles without Andreev re- 
flections and show how it is related to the density of 
states of the superconductor. We will then incorporate 
Andreev reflection using the BTK model, and show how 
to estimate the energy gap and place a lower bound on 
the Fermi velocity and BCS coherence length. 

Ignoring Andreev reflections for the moment, the prob- 
ability for a tunneling process is the product of the num- 
ber of occupied initial states, the number of unoccupied 
final states, and the probability that the tunneling event 
occurs. For tunneling between the normal metal to the 
superconductor, the number of occupied initial states is 
Nn(E)f(E) and the number of unoccupied final states is 
N s (E + eV)[l- f(E+eV)]. Here, N N {E) is the density of 
states of the normal metal, Nq(E) is the density of states 
of the superconductor, and f(E) is the Fermi function. 
For the reverse current, the joint probability for this pro- 
cess is N N (E)f(E + eV)N s (E + eV)[l - f{E)]. Hence if 
we subtract the reverse current from the forward current, 
we obtain 



Ins = a|T| 2 / N N (E)N s (E+eV)[f(E)-f(E+eV)} dE, 

J — OO 

(20) 

where a is a constant that depends on the area of the 
junction and the Fermi velocity (see pp. 44-47 of Ref. 
12). \T\ 2 is a phenomenological tunneling matrix ele- 
ment that depends on the type of insulator between the 
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N/S junction and governs whether the event occurs. The 
difference in Fermi functions f(E) — f(E+eV) represents 
the probability of transmission, where V is the applied 
voltage and eV the resulting difference in the chemical 
potential across the junction. 

In the energy regime in which we are dealing, we can 
approximate N-^(E) = N^(0) as constant since E <C M 
has been previously assumed. Using the current in Eq. 
(20), the differential conductance is then 
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(21) 



where we have assumed that N§(E) is slowly varying 
compared to the Fermi weighting factor. Note that 
—df(E + eV)/d(eV) is a bell-shaped weighting func- 
tion with a maxinum at E = — eU, a width of ~ k^T, 
and unit area. This factor plays the same mathemati- 
cal role as a Dirac delta function inside an integral when 
k^T — > 0. At zero temperature, the differential conduc- 
tance approaches the density of states of the supercon- 
ductor 
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which is equivalent to what is shown in Fig. 2b for large 
Z. 

If we use the BTK formulation to incorporate Andreev 
reflections, the current that flows across the junction is 

/oo 
(1 + A(E) - B(E))(f(E) - f(E + eV))dE. 
-OO 

(23) 

In this equation, 1 + A(E) — B(E) represents that trans- 
mission coefficient for electrical current. This coefficient 
incorporates the fact that when a hole is reflected, a posi- 
tive charge moves in the opposite direction as the incident 
electron. Such a current is represented by A{E) and will 
add to the incident current. Similarly, reflected electrons 
are represented by B(E) and will reduce the current. Al- 
though the current obtained in Eq. (23) using the BTK 
approach is similar to Eq. (20), it is different in the fact 
that the normal metal and superconductor are not de- 
coupled and the initial wave packet evolves continuously 
through the barrier (see section V of Ref. 3). However, 
the first approach is important from a pedagogical stand- 
point, as it is easy to interpret the relationship between 
the tunneling current and the density of states of the 
superconductor. 

Using the current in Eq. (23), the differential conduc- 
tance of a S/N junction at low temperature is approxi- 
mately given by the transmission coefficient for electrical 
current 
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d^NS 



NS 



dV 



oc (l + A(eV) -B(eV)). 



(24) 



T=0 



This conductance is also related to the density of states 
of the superconductor, but the relationship is harder to 



discern because the density of states is encoded in the 
coherence factors N$(E) — (uq — i>o) _1 - A benefit to 
this BTK formulation is that the barrier height is not 
assumed to be large so that Andreev reflection processes 
come into consideration, whereas it is absent in the for- 
mulation given by the current in Eq. (20). Therefore, 
Eq. (24) provides a spectrum of conductances from a low- 
Z barrier (Fig. 2a) to a high-Z barrier (Fig. 2b). 

By measuring the differential conductance of a Nb/Au 
junction and fitting the experimental data to the BTK 
model using the derivative of Eq. (23), we can obtain val- 
ues for Z c ff , A, and T. The latter two are of fundamental 
interest. In addition, we can use the value for Z e g to 
obtain a lower bound on the Fermi velocity in the super- 
conductor. Letting Z — and using Eq. (19) to obtain 
a lower bound on the Fermi velocity ratio r, we find the 
BCS coherence length £ of the quasiparticle waves in Nb 6 
to be 

Hvf s 



7rA ' 



(25) 



It is important to note that the BTK model assumes 
ballistic transport, which implies that the size of the con- 
tact a is less than the electron mean free path I and the 
coherence length £. If a <C I, the contact is ballistic 
which prevents heating effects at large current densities 
at voltages on the order of the gap energy. The resis- 
tance in this regime, known as the Sharvin resistance, 
Rq = pi/ 4:0? , is due to the radius of the constriction from 
the normal metal to the superconductor. Furthermore, 
the size of the contact must be much less than the co- 
herence length (a <C £) to suppress the proximity effect, 
where the tip of the gold begins to weakly superconduct 
and the niobium becomes an imperfect superconductor. 
Additionally, this condition prevents the destruction of 
superconductivity near the gap where the carrier veloc- 
ity approaches the depairing value (A/pp), where pf is 
the Fermi momentum. 6 

In a typical low temperature superconductor, the co- 
herence length is on the order of 1000 A and the contact 
in the Sharvin regime is around a ~ 100 A. However, typ- 
ical tips, such as a gold wire cut by a razor blade, or in 
our case an electrochemically etched gold tip, have a ra- 
dius on the order of 0.1 — 1 /zm. A tip this blunt should 
have a contact resistance around 1 — 10 mil. In prac- 
tice, the contact resistance is often found to be around 
10 — 100 il, corresponding to a much smaller contact area. 
Various explanations have been proposed for this obser- 
vation, but there is no consensus as to why the effective 
contact size is so good in practice. Nevertheless, data 
can be found that fits well with BTK theory. 6 



C. Fitting Experimental Data to the BTK Model 

We now explain how to fit the theoretical model to ob- 
tain the parameters A, T, and Z c g . In Eq. (24), the con- 
stant of proportionality is unknown, so in order to com- 
pare the experimental and theoretical results, the data 
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is normalized by a value of the differential conductance 
above the gap. It is not uncommon for there to exist a 
background slope due to thermal voltages or other exter- 
nal factors. To try to minimize this effect, we applied the 
following normalization convention 
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which is shown as a dashed line in Fig. 2. 

In order to perform a non-linear least squares fit of 
the data to the theoretical model, an array of theoretical 
conductances given measured voltages needs to be con- 
structed. In our case, the temperatures are significantly 
below T c , so we can approximate the differential conduc- 
tance as in Eq. (24). It is convenient to use MATLAB 
or some other program to expand the quantities in this 
equation, particularly since a and b are complex. The 
measured voltage values are then taken as the input for 
the theoretical differential conductance Gthcor(eV) in Eq. 
(24). 

For cases where the temperature is closer to T c , we 
can perform a Riemann sum and cut off the limits of 
integration because the integrand is non-zero for a small 
range localized about eV. Both methods were used for 
our data analysis in order to confirm that allowing T = 
did not alter the values of the parameters. 

Together, the values of A, T, and Z c g parameterize the 
differential conductance of Eq. (24). An additional pa- 
rameter was included to account for any voltage offest. A 
nonlinear least squares curve fit with MATLAB was used 
to find the values of these parameters where the theoret- 
ical and experimental curves most closely matched. 



III. EXPERIMENTAL SETUP 

The Nb/Au junction is formed by having a differen- 
tial screw that allows an electrochemically etched Au tip 
to gradually touch down onto a Nb film. The Nb film 
was made in an evaporator that was baked overnight to 
minimize the amount of water in the chamber, as excess 
oxygen during deposition decreases the transition tem- 
perature. A silicon substrate was ion milled for a few 
seconds in order to remove excess water or contaminants. 
Then 100 nm of Nb was evaporated onto the substrate 
at a rate of 2 A/s at a base pressure of ~ 10~ 7 Torr. We 
independently verified a transition temperature T c ps 7 
K for the sample using a four-point measurement of re- 
sistance vs. temperature. 

The sample was adhered with Lakeshore VGE-7031 in- 
sulating varnish to a chip holder made from a small PC 
board affixed with Loctite Hysol 1C epoxy to a copper 
beryllium plate. Wire bonds were used to make electrical 
contact between the PC board and the sample. To mea- 
sure the differential conductance through the junction, 
two wires from the sample are necessary for a voltage 
and a current lead. However, four contacts were placed 
on the sample to perform a four point measurement of 
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FIG. 5: (Color online) A cross-section of the bottom end of 
the cylindrical probe (not to scale). All disks, rods, and holes 
are cylindrical and all parts are fixed except for the inner rod 
(dark gray). The outer rod (light gray) has a smaller diameter 
at the end in order to reduce horizontal motion of the inner 
rod. The Nb sample is on a removable plate that is varnished 
to the bottom disk. Wires from both the Nb and Au are 
threaded up the probe through a hole cut out in the upper 
ring; each wire drawn corresponds to a pair of wires threaded 
around each other. 



the resistance to verify that the sample became supercon- 
ducting when placed in liquid helium. We then varnished 
the copper-beryllium plate to the base of our apparatus. 

The gold tip was electrochemically etched 18-22 start- 
ing from a 0.5 mm gold wire of 99.9985% purity pur- 
chased from Alpha Aesar (product number 10966). The 
tip acted as an anode while a platinum foil acted as a 
cathode in a diluted 40% HC1 solution that was heated 
to 50°C. A function generator produced a 10 kHz square- 
wave signal with 5 V peak-to-peak and 2.5 V offset in 
order to act as the voltage source for the electrochemi- 
cal reaction. Immediately after etching, we rinsed the tip 
with hot deionized water to remove salts and then dipped 
it into acetone for further cleaning. We then placed the 
gold tip into an insulating acrylic plastic holder attached 
to the apparatus (Fig. 5). While the plastic holder func- 
tioned well for our purposes, it did crack slightly. For 
future work we recommend using macor or polycarbon- 
ate, both of which are insulating machinable materials 
that have better thermal contraction properties at cryo- 
genic temperatures. 

A model DM-13 differential micrometer purchased 
from Newport controlled the vertical movement of the 
tip. It was mounted at the top of the probe outside the 
glass dewar cryostat and on top of springs that mitigated 
vibrations (Fig. 6). An O-ring seal held the dipping 
probe in place, but the force supplied by the differential 
screw was sufficient to overcome this grip and allow us to 
push the probe further toward the sample. The O-ring 
also provided a good enough vacuum to prevent the es- 
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FIG. 6: (Color online) A cross-section of the top part of the 
probe (not to scale). All holes and disks are cylindrical, with 
radii and heights specified where appropriate. The inner most 
component of the micrometer is attached to the inner rod 
(light gray), which is free to move vertically. The inner rod 
is attached to the central plate which is also free to move. 
Springs attached to the three support rods help to dampen 
noise. A nut with an O-ring seal, used in ultra-torr vacuum 
fittings, is used to loosely hold the inner rod in place with re- 
spect to the outer rod (dark gray) and aids in holding vacuum 
while measurements are taken. The wires are epoxied into a 
removable brass tube which is then inserted into a permanent 
tube and sealed with a nut and an O-ring. When pumping on 
the liquid helium, a rubber stopper is inserted into the hole 
used for transferring helium. The bottom plate is screwed 
into the top of the dewar. 



cape of helium gas from the bath. In order to minimize 
unwanted horizontal movement of the tip, we welded an 
outer rod to the top plate to act as a fitted guide for the 
inner rod driven by the differential screw. There was a 
guided bearing at the base of the outer rod that mini- 
mized the horizontal motion of the inner rod (see Fig. 
5). Stainless steel was used to construct the inner and 
outer rods. See supplementary material for more detail 
on the schematics of this apparatus and the cost of the 
electronics. 23 

The differential conductance was obtained for differ- 
ent input currents by using an AC modulating technique 
(p. 46 of Ref. 24). A DC current 7 D c and an AC cur- 
rent Jac cos(wi) with Jac "C Ibc was applied across the 
junction and the voltage was measured. Expanding the 
voltage in a Taylor series 



^(/dc+/ac cos(wi)) 
dV(Ipc) 
dl 



V(Iuc) 

I AC COs(ujt) 



(26) 



we see that a measurement of the AC voltage determines 
dV/dl, the inverse of the differential conductance. 

Figure 7 depicts a block diagram of the electronics. A 
Keithley 6221 AC/DC current source provided the AC 
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amplifier 



current source 



lock-in 



L| PC board ^ 
LHe 

dewar 



multimeter 



FIG. 7: (Color online) A block diagram of the electronics. 



and DC currents across the junction. We used a lock-in 
amplifier across the junction to measure the AC voltage 
and a digital multimeter to measure the DC voltage. For 
contact resistances from 2-20 f2, an AC current of around 
50 fj,A produced an AC voltage approximately an order 
of magnitude less than the gap for Nb (~ 1 meV). In 
general the AC voltage used should be about an order 
of magnitude less than the width of the peak at the gap 
energy in order to resolve the signal. If the AC voltage 
were larger than the peak width, then the differential 
conductance at the gap energy would be in essence an 
average of measurements not localized around the peak 
height, leading to an underestimate of its true value. We 
used a frequency of 1 kHz for the AC current and the 
DC current was stepped from negative to positive values 
and back so as to observe any heating effects or other 
anomalies. For most measurements |/dc| < 600 /iA. 

When the gold tip approaches the sample, a resistor in 
parallel with the source will prevent overloading the cur- 
rent source. This resistor should be a few orders of mag- 
nitude greater than the expected resistance of the N/S 
junction. Once the tip touches down onto the sample, 
the resistance will drop significantly, indicating that con- 
tact is made. At this point, the shunt resistor is removed 
from the circuit. In addition to the electrical signal that 
is observed when contact is made, the glass dewar also 
has a slit to visually observe when the tip comes in con- 
tact with the sample. The differential screw allows us 
to increase and decrease the pressure of the tip on the 
sample, thus changing the contact resistance for various 
measurements. 

While it is possible to have the Au tip in contact with 
the Nb film at room temperature, this can cause thermal 
contractions that damage the sharp tip and increase the 
contact area beyond that required to observe Andrccv 
reflection. Therefore, the tip and sample are separated 
at the start of the experiment. 
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FIG. 8: (Color online) Normalized conductance with dots rep- 
resenting experimental data and solid curves representing a fit 
to BTK theory, (a) A sweep from negative to positive bias 
voltage with VoSset = 0.33 mV at 10-f2 contact resistance, (b) 
A sweep from positive to negative bias voltage at a higher tip 
pressure, with T4ffset = —0.11 mV at 3-fi contact resistance. 





io n 


3 n 


Park et. al. 


A (meV) 
r (meV) 


0.82 ±0.04 
1.53 ±0.10 
0.78 ±0.10 


0.29 ±0.03 
0.80 ±0.01 
0.02 ±0.01 


0.453 
1.25 
0.55 



TABLE I: Comparison of Z e s, A and F as determined exper- 
imentally with both a 10 fi and 3 Q contact resistance. 



IV. RESULTS 

Most of the our data is similar in character to the mea- 
surement shown in Fig. 8a (for contact resistances around 
10 f2) and Fig. 8b (for contact resistances around 3 Q). 
These measurements were performed at a temperature of 
T ~ 1.5 K. Table I summarizes the fit parameters Z e ff, A 
and T, which were found using a MATLAB least squares 
fit to Eq. (24). For comparison, the fit parameters from 
Park et. al. are also provided. 5 

For the data shown in Fig. 8a, the current was swept 
from negative to positive bias voltage and a DC offset of 
0.33 mV was used to shift the original data to be cen- 





io n 


Park et. al. 


Tc (K) 
v Fs (m/s) 
£ (nm) 


7 

> 3.1 x 10 5 
43 


9.22 
> 5.8 x 10 5 
97 



TABLE II: Comparison of estimated values of vf s and £ for 
Nb. 



tered around zero voltage. The offset used in fitting data 
for the sweep back from positive to negative bias was 
w —0.11 mV. The reason for this offset is unknown, but 
could be due to an offset in our electronics or the presence 
of a time constant. From Z e ff , we note that the differen- 
tial conductance is neither fully in the Andreev reflection 
regime nor fully in the tunneling regime, but instead is 
a combination of the two. This measurement was made 
after decreasing the tip pressure to increase the contact 
resistance. 

Given our effective barrier strength Z c g and the Fermi 
velocity of Au («f n = l-4x 10 6 m/s), 6 we estimate a lower 
bound on the Fermi velocity of Nb to be 3.1 x 10 5 m/s 
and a BCS coherence length of 43 nm. For comparison, 
an experiment that measured the superconducting upper 
critical field H C 2 in bulk polycrystals found the Fermi 
velocity of Nb to be vp s — 5.7 x 10 5 m/s, consistent with 
our results. 25 Additionally, Park et. al. performed point 
contact spectroscopy with a gold tip on a 210-nm film of 
Nb with T c — 9.22 K and obtained similar results after 
fitting their data to BTK theory. 5 The pertinent results 
are listed in Tables I and II for comparison. 

As shown in Table II, our measured coherence length 
is noticeably smaller than that reported by Park et. al. 
Because of the impurities in our Nb thin film that reduce 
T c to 7 K, we expect a smaller coherence length than for 
a pure Nb film with T c of 9.2 K. Other experiments have 
found the coherence length for bulk Nb to be around 40 
nm. 26 Hence our measurements are in reasonable agree- 
ment with the literature. 

After pressing the tip further into the sample to reduce 
the contact resistance to 3 we obtain significantly re- 
duced values for A and T (see Fig. 8b and Table I) and 
a dip in the differential conductance at around ±2 meV. 
One possible explanation for the reduced gap is the prox- 
imity effect, in which A(x) departs from simple step func- 
tion behavior at the interface of S/N, and the end of the 
probe tip becomes weakly superconductivity. The prox- 
imity effect is known to decrease the observed value for 
A. 5,27-29 Furthermore, the dip in the conductance was 
seen previously in a Nb/Cu point-contact spectroscopy 
experiment by Blonder and Tinkham who found that the 
tip had curled back on itself. Blonder and Tinkham ar- 
gue that this wider contact area allows heating to play a 
significant role. 16 At this point, the initial assumptions 
of the BTK model — that a <C £, £ and that the electrons 
are ballistic — begin to break down. 6 ' 30 We also observed 
a curled tip in our experiment. 
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V. CONCLUSION 

We performed point-contact spectroscopy to measure 
the differential conductance as a function of bias- voltage 
through a Au/Nb junction. We then estimated the su- 
perconducting gap, quasiparticle lifetime, and placed a 
lower bound on the Fermi velocity and BCS coherence 
length in Nb by fitting the experimental data to theoret- 
ical curves predicted by the BTK model. 

One of the largest drawbacks to this setup is the 
amount of time required to cool down a glass dewar cryo- 
stat before data can be acquired. In addition, we only ob- 
tained approximately ten reasonable sets of data before 
a dip in the differential conductance appeared around 
±2 mV. As discussed above, this signifies a curled tip 
with large contact area that negates the underlying as- 
sumptions of the BTK model. Unfortunately, at this 
point a new tip cannot be easily replaced mid-experiment 
because the cryostat must be warmed up to remove the 
tip. Furthermore, there were concerns about the possi- 
bility of the probe freezing due to thermal contraction of 
the cylindrical rods. Such contraction would trap helium 
between the rods and render the probe immovable. Al- 
though such concerns were never realized in practice, this 
possibility could be eliminated by cutting a small slit in 
the outer rod to allow helium to escape. 

This experiment easily extends to measure other super- 
conducting materials or different types of junctions, such 
as a superconductor-superconductor junction. Plus, the 
introduction of a magnetic field into the apparatus would 
allow for measurement of tunneling processes from a nor- 



mal metal to a ferromagnet. One improvement would be 
to measure incremental changes on the differential screw 
to provide quantitative data of the pressure between the 
tip and the sample. 

The present experiment can be simplified by using a 
block of Nb rather than a thin film. Films are easier to 
insulate and were readily available to us, but a Nb block 
can be placed on a glass slide and wires attached with 
silver epoxy. A superconducting sample can also be pur- 
chased from industry. If one obtained a high-T c sample 
(e.g. YBa2Cu307_5), a simple setup can be constructed 
to immerse the sample into a liquid nitrogen container. If 
liquid helium is required, one can reconfigure the appara- 
tus described previously to fit inside of a storage dewar. 
Because the height at which the dipping probe is inserted 
into the storage dewar stabilizes to a fixed temperature, 
such a setup would allow measurement of the differen- 
tial conductance as a function of temperature. Hence 
the transition of the differential conductance when the 
sample is above T c to below T c can be observed. 
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